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^ ; Abstract. 

^T) • Equation for the sum of BFKL pomeron fan diagrams is rederived by direct summation 

^ , and solved numerically for rapidities y < 50. At high rapidities y > 20 the resulting cross- 

sections for the scattering of a longitudinally polarized qq pair on the nucleus cease to depend 
' on its transverse dimension and tend to a constant limit 0.1768 R\, which corresponds to 

scattering of a colour dipole on a black disk. Thus the unitarity is restored and the singularity 
in the j plane is reduced to a simple pole at j = 1. The nuclear structure function at small 
■ X behaves as ln(l/x). The found gluon density has a soliton-like form in the log/c space: 

. its form is close to Gaussian, independent of rapidity, the center moving towards higher log k 

Qh! with a nearly constant velocity as rapidity increases. 

Oh: 
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In the framework of the colour dipole model of A.H.Mueller [1,2] it follows that in the high- 
colour limit Nc ^ oo the scattering on a heavy nucleus is exactly described by the sum of 
' fan diagrams constructed of BFKL pomerons, each of them splitting into two [3]. This sum 

seems to be unitary by itself. It is important that no splitting into three or more pomerons, 
as introduced in [4], occurs, although formally they give contributions of the same order. 
Because of this fact, once the splitting vertex is known, construction of the amplitude for the 
interaction with the nucleus becomes straightforward, reducing to summing BFKL pomeron 
fan diagrams. This procedure has been well-known since the times of the old Regge-Gribov 
theory [5]. 

In the perurbative QCD a pioneering step was taken also many years ago by L.Gribov, 
E. Levin and M.Ryskin, who summed fan diagrams in the double log approximation and wrote 
their well-known non-linear GLR equation [6]. In the framework of the BFKL dynamics the 
necessary tool for constructing fan diagrams is the corresponding triple pomeron vertex, which 
was found in the colour diploic approach for Nc ^ oo hy A.H.Mueller and B.Patel [2] and 
in the s-channel unitarity approach for any number of colours by J.Bartels and M.Wuesthoff 
[7]. The equivalence of both results was shown in [8]. 

The equation for the sum of BFKL fan diagrams with this splitting vertex was written 
by I.Balitsky [9] in his original operator expansion formalism and by Yu.Kovchegov [10] in 
the colour dipole framework (with a somewhat unconventional form of the coupling to the 
target) . Its perturbative solution in the region of small non-linearity (outside the saturation 
region) was studied in [11]. Asymptotic estimates of the solution were presented in [12]. 
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In the present paper we first rederive the BFKL fan diagram equation by direct summation 
using the standard form of the pomeron-target couphng. It has a form of a simple (and 
elegant) evolution equation in rapidity y for a wave function ^{y, q) in the momentum space: 

'-^ = -H4>{y,q)-4>\y,,l (1) 

where H is the BFKL Hamiltonian for the so-called semi-amputated function (a similar form 
was also obtained in [11]). 

In spite of its tantalizing simplicity, Eq. (1) does not seem to allow for an analytical 
treatment except by perturbative methods, not valid in the most interesting region of strong 
non-linear effects, or by qualitative asymptotic estimates. The bulk of this paper is corre- 
spondingly devoted to its numerical analysis. We numerically study evolution of the wave 
functuion in rapidity, starting from an appropriately chosen initial function. The results are 
then used to find the structure function of the nucleus at small x = and various vir- 
tualities Q^. Our results show that for a heavy nucleus the longitudinal part of structure 
function saturates at x — > to a universal function F2^\q'^), independent of the nucleus 
atomic number but strongly dependent on Q"^ in the whole range of < 10^ (GeV/c)^ 
explored. In fact pi^^Q"^) is just proportional to Q^. This reflects the behaviour of the 
scattering cross -section of a longitudinally polarized qq pair on the nucleus: at very small x 
it becomes independent of both the pair size and x. The latter property indicates that the 
unitarity is restored and the leading singularity in the complex momentum j is reduced to a 
simple pole at j = 1. Note that the limiting cross-section is in accordance with a picture in 
which the longitudinally polarized photon with a certain probability splits into colour dipoles, 
which then scatter on the nucleus as on a black disk. For the leading transverse part of the 
structure function this probability results infinite. Due to this fact both the cross-section and 
the structure function continue to grow at a; ^ approximately as ln(l/a;). These results 
are in agreement with predictions made in [11,12]. 

As to the gluon density, we have found it to have a Gaussian shape as a function of 
^ = In A; with a center at y. So with the growth of rapidity it propagates towards higher 
momenta, practically preserving its form, very much like a soliton wave. Obvious limitations 
on computing time and memory allow to follow its movement up to momenta not higher than 
of the order lO^'^ GeV/c. However we expect that this behaviour persists in the model until 
arbitrary high values of momenta. As a result, at any fixed value of k the density eventually 
goes to zero as y ^ oo. In this sense we have the maximal saturation of the density possible. 

The paper is organized as follows. In section 2 we present derivation of the BFKL fan 
diagram equation by direct summation. Section 3 is devoted to the numerical solution of Eq. 
(1). In Section 4 the nuclear structure function and gluon density are calculated. Section 5 
contains our conclusions. In the Appendix we compare our BFKL fan diagram equation with 
the Kovchegov's one [10]. 



2 Fan diagram equation for the BFKL pomerons 

We start with a single scattering contribution to the forward ampitudc for the interaction of 
the projectile particle with a nucleus. In the BFKL framework, at fixed impact parameter 6, 
it has a well-known form 



Ai{y,h) = isg^AT{h) / d\d\' p{r)G{y,r,r')pN{r'). 



(2) 
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Here p and pN are the colour densities of the projectile and the target nucleon respectively. 
Function G is the forwrad BFKL Green function [13]: 

where and (p' are the azimuthal angles and 

uj{u) = 2(asiV/27r)(V'(l) - Re^(l/2 + iu)) (4) 

are the BFKL levels. Due to the azimuthal symmetry of the projectile colour density one 
may retain only the term with zero orbital momenta n = in (3). Separating the projectile 
part, the single scattering term may be written in the form 

Ai{y,b) = 2is J d^rp{r)^^{y,b,r), (5) 

where 

b, r) = ^g^AT{b) J d\'G{y, r, r')piv(r'). (6) 

The double scattering contribution has been calculated in [14]. In the limit Nc» 1 one 
finds 

A2{y,b) = -'^^-^ J d^rpir) dy2 J '[[d^riS'^iri + r2 + r^) 
r?Vf 

-^-^G{y - y2, r, ri)$i(y2, b, r2)$i(y2, b, -rs). (7) 



^2^3 

Presenting A2 in terms of $2 similarly to (5) we find 

8^ 



$2(y, b, r) = / dy2 11 d\d\ri + + rg) 

^0 J -^^ 



r?Vf 

-^r^G{y - y2, r, ri)$i(7/2, b, r2)$i(y2, b, -rs). (8) 



^2^3 



The whole set of fan diagrams will be evidently summed by the equation which is graph- 
ically illustrated in Fig. 1: 

$(y, 6, r) = $i(y, b, r) - j dy2 J ]]_ d\S^{ri + r2 + r^) 

^ 1=1 

-\^G{y - y2, r, ri)$(j/2, b, r2)$(y2, 6, -^3). (9) 



The impact parameter b appears here only as a parameter and the dependence on it will be 
implicit in the following. In terms of $ the total forward scattering amplitude on the nucleus 
will be given at fixed b by 

Aiy,b) = 2is j d^rp{r)^{y,b,r). (10) 

One can rewrite Eq. (9) as an evolution equation in y. To this end we present $ as an 
integral over i/'s similar to (3): 



$(y,r) = j dvr^-'^'''^{y,v). 
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Then from (9) we find 

.4 



= 64.^(J+1/4)^ ^^^^^V dW^-'^-'Pir'^^nb)- 

32.5(:.;i/4)^ °"""' /. / n ""^--.^^in + --^ + '-3) 



'■i5r!-«"*(!,2,r2)4.te, -,-3). (11) 



Tor 
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Multiplying both parts by e '^'^^^^ and taking a derivative in y we obtain 



^r|-2-$(y,r2)$(y,-r3). (12) 
Returning to the r space we take into account that 

u{iyy-^''' = -Hr^-'^'\ (13) 
where H is the BFKL Hamiltonian [13]. Then (12) transforms into 

— + HU{y,r) = -^ / ]ld\S''{n+r2+rs)^G{0,r,nMy,r2my,-rs). (14) 



i^— + Hj^y,r) = -^ J '[ld\d\r + r2 + rs)-j-i^{y,r2My,-r3). (16) 



Using 

ll^G(0,r,ri) = -^ ^ e-('^-'^i) / dKr/ri)-^^^ = <52(r - n), (15) 
we simplify (14) to a form 

^ (fri5^{r + r2 + r-a)-^-? 

=2 ^2^3 

The initial condition is determined from (9) to be 

^(y,Oy=o = *o(r) (17) 

with 

$o(r) = \g^AT{h) I d\'G{0,r,r')pir'). (18) 

Except for the initial condition and different variables, Eq. (16) coincides with the one 
constructed by Kovchegov in the colour dipole approach in [10] (see Appendix). 
Now we go to function 

0(y,r) = l$(y,r) (19) 
and pass to the momentum space. For <^{y, q) Eq. (16) reads 
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where we have used that (f){y,q) in fact depends only on \q\, which follows from the initial 
conditon 



My,(l)y=o = J ^e-"^'^o{r) (21) 

with "3>o given by (18) and depending only on |r| due to the azimuthal summetry. 

The Hamiltonial r^'^Hr'^ which appears on the left-hand side is the standard forward 
Hamiltonian for the semi-amputated functions 

liJr^ = H = + lnr-2 - 2(ln2 + V(l))]. (22) 

Indeed the eigenfunctions of H = r^-'^'^^ e^'^^ go over into the eigenfunctions of H 

after division by r^. 

This brings us to nearly the final form of our equation 

f I- + h] 0(y, q) = -4^<PHy, q). (23) 



It remains only to appropriately rescale H, (j) and y to obtain the BFKL fan diagram equation 
in the final form (1). We introduce 

H = ^H, y=^y, and cP{y,q) = 27T4>{y,q), (24) 

where 

H = lnr^ + lnq'^ - 2(ln 2 + '0(1)). (25) 
The equation for <f) takes the form (1) 

(^^ + Ho^Hy,Q) = -$\y,Q) (26) 

with the initial condition 

4>iy,Q)y=0 = ^^^^ I dPrcPn^e^'i^G{^,r,n)pN{n). (27) 

3 Numerical solution of the BFKL fan diagram equation 

To solve numerically Eq. (26) wc first transform the BFKL Hamiltonian H to more convenient 
variables. In the momentum space one can write action of on a function ^(gf) as 



oo 



H<j^iq) = -2l kdk ' .J ^^^^/--v.. ^^ j..,^ ^ _ (28) 



We transform it to variables u,v which take values in [0,1]: 

(J = cxp[(Mi +M2)w - Ml], fc = exp[(Mi +M2)u - Ml], (29) 

where q, k are in GeV/c and Mi(2) is a lower(upper) integration limit in In A;. Then (28) goes 
into 

"1 r/A(^) _ J(^ _ ^)(^(|i) i^(^t) 



H<p{u) = -2(Mi + M2) / dv ' . ' + . , ' . , (so) 

h L \f(u-v)-\\ + ^j-\u-vy 

where 

/(«) = exp[(Mi +M2)u]. (31) 
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Now we discretize the interval [0,1] in u and v into n equidistant points uo,ui, ...Un and 
vq, vi, ...Vn to convert the evolution equation (26) into a set of n ordinary 1st order non-linear 
differential equations in y. This set of equations can be solved by standard methods. We used 
the simplest 2nd order Runge-Kutta algorithm. The standard values of n and of the number 
of iterations to evolve in two units of y were 800. We have studied evolution of (p{y, q) from 
the initial value at y = up to y = 10. The (fixed) value of the strong coupling ag = 3^/47r 
has been chosen to be 0.2. With this choice our maximal rapidity is around 50. 

An evident difficulty which one meets are the values of H(f) at endpoints uq and n„, at 
which the introduced cutolfs make the results not reliable. To overcome this difficulty we 
calculated Hcf) at these points by extrapolation from the neighbouring points. The stability 
of this procedure was checked by comparing the results for double and quadruple values of n. 

The initial function is determined by the colour density of the nucleon pN{f) according 
to (27). To simplify our calculations we have taken the Yukawa form for pn{t): 

where //=l/0.7 fm has a meaning of the inverse nucleon radius. This choice may look a bit 
arbitrary, but our results show that with the growing rapidity the system quickly forgets 
not only the form but even the absolute magnitude of pM-, so that the solution becomes 
independent of the initial function and governed only by the internal dynamics of Eq. (26) 
itself. 

Doing the integrations over r and ri in (27) and over v inside the BFKL Green function 
(3) we find ioi q > p, 



and ioT q < p 

Mq)<i<^ = b{\[{2 In ^ + V'(3/2) - ^(1) - 1)' + ^-'(3/2) + ^-'(1) + 1]- 



n=0 

where the dimensionless coefficient 



.4- p 

_^V(-l)"f — 1 "^"^"^"^ \ (34) 




2 (35) 



167r p, 

carries all the information about the nucleus. Its maximal value with ag = 0.2 is about 0.12 
for the central scattering on lead (b = 0). The asymptotic behaviour of (j)o{q) at (7 ^ oo and 
g ^ is governed by the first terms in (33) and (34), which come from the poles of the BFKL 
Green functions at v = zbz/2: 

^o(g) ~ g^oo, 

(^o(g) -log'?, q^O. (36) 

We have checked that our results parctically do not change if the complicated function which 
multiplies B in (33) and (38) is substituted just by ln^(g/(/x + q) with the same asymptotic 
behaviour (36). 

The initial function depends on the form of the nucleus profile function T(b). For our 
numerical calculation we have chosen 



T{b) = 2JR\ - h^/VA (37) 
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corresponding to a finite nucleus of radius Ra and volume Va with a constant density. We 
have taken Ra = A'-^^Rq with Ro = 1.2 fm. 

Our results for (j){q) for two values of i? = 0.12 and 0.02, corresponding to a central and 
very peripheral collisions off lead, respectively, arc presented in Figs 2,3. One observes that 
the solution cp at rapidities y > 2 evolves to a very simple and universal form, practically 
independent of the initial function. Crudely speaking it linearly falls with ^ = Ink until it 
meets the x-axis wherefrom it stays equal to zero. The slope of the falling part is exactly 
equal to 1, so that very crudely 

m = Uy,b)-^, for e<Ci(y,^'), m = o, fore>a(y,&)- (38) 

The value of (y, 6) and hence the interception point with the x axis grow linearly with y, 
so that with the growth of y the picture simply shifts to the right. In reality the curve for 

4> of course has no break: the two straight lines of which it is formed join smoothly in the 
vicinity of ^i. As will be clear later, the physically important region is precisely this vicinity, 
where (p is not trivial. 

The results for <^{k) do not depend on the chosen cutoffs, provided they are taken to 
cover the region around the interception point ^i. Otherwise the evolution stops as soon as 
^1 touches the upper cutoff M2. So if one wants to study evolution up to high values of y the 
upper cutoff should be taken correspondingly high. In our calculations we chose Mi = 10 
and M2 = 20, having verified that further raising of either Mi or M2 does not change the 
results. 

In conclusion we find that (f) docs not possess any finite limit as y ^ 00 so that Eq. (26) 
does not lead to any saturation of the wave function at high rapidities, contrary to naive 
expectations. However in the next chapter we shall see, that such saturation indeed occurs 
for physical quantities. The point is that function (j) by itself has no physical meaning. It is 
its derivatives which matter. 



4 Nuclear structure function and the gluon density 

The nuclear structure function is obtained in the standard manner as 

F2{x,Q^) = ^{aT + aL), (39) 

where (Tt,l arc the total cross-section for the scattering on the nucleus of a virtual photon 
with transversal (T) or longidunal (L) polarization. Both cross-sections can be found from 
the imaginary part of the forward scattering amplitude (10). In terms of (p we find 

aT,L = y (fb(frpT,L{r)r^Hy, r, b). (40) 

Here we explicitly indicated the dependence of cp on the impact parameter; pT,L{r) are the 
well-known colour densities of the virtual photon split into a qq pair (see e.g. [15]). With 
massless quarks 

PTir) = da{a' + (1 - af)e'Kl{er)) (41) 

and 

Pl{t) = ^^Q' I daa\l - afKl{er), (42) 
where = Q'^a{l — a) and is a sum of squares of quark electric charges in units e. 
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Passing to momentum space we find 

/(Pq ~ 
d^bj^;^4>{q,y,b)wT,L{q), (43) 

where 

wt,l{q) = I d\r^pT,L{r)e"ir (44) 

Straightforward calculation leads to the following expressions for wt,l{<i)- For the transverse 
density one findg|^ 

^T{q) = ^^ da{a^ + {l-af)Vl[{q^/2 + e^)J{q,e)], (45) 

where 



The longitudinal density w/, is given by the same expression with substitutions 

+ {l-af ^ a{l-a), q^/2 + e^ ^ -Ae^ 

With the found numerical values for the function </> in the range < y < 10 we evaluated 
the nuclear structure function of lead {A = 207) for various 

Q2 between 3 and 10^ (GeV/c)2. 
The results are most instructive for the the cross-section gl for the scattering of a longitu- 
dinally polarized virtual photon on the nucleus. In Figs. 4,5 we present it (with 1) 
as a function of y at fixed and as a function of at fixed y respectively (in GeV~^). 
From Fig. 4 one clearly sees saturation in rapidity: for any value of the cross-section 
tends to the same limit of 228.8 (GeV/c)~^ = 0.1768i?^ as y goes beyond 5. The resulting 
constant cross-section is evidently consistent with the unitarity restrictions. In terms of the 
complex angular momentum j our results indicate that the original cut at j > 1 is reduced 
to a simple pole at j = 1. Fig. 5 illustrates an unusual behaviour in which sets in at high 
y: instead of going down as \/Q in the standard BFKL approach, the cross-section becomes 
independent of to a very high precision. 

These results for the longitudinal cross-section can be conveniently interpreted in terms of 
scattering of colour dipoles off the nucleus. The density piif), appropriately normalized, can 
be interpreted as a probability distribution for the longitudinal photon to split into colour 
dipoles of transverse dimension r. The normalization factor 



r p2 TU z2 

D= d^rpiir) = -j^^ = 0.028145 



can be considered as a total probability for the photon to split into dipoles. Then the dipole- 
nucleus cross-section is found by dividing by factor D, which gives 8131.1 (GeV/c)~^ 
independent of Q^, that is, of the dipole dimension. This value exactly equals 2'kR\, corre- 
sponding to scattering off a black disk. 

The transverse part of the structure function does not admit this interpretation, due to 
the fact that PT{f) is not normalizable. As a result both the cross-section and the structure 
function do not saturate at large y but continue to grow nearly linearly in y. This is illustrated 
in Figs. 6,7 where the total structure function (including the much smaller longitudinal part) 
is shown as a function of a; = and respectively. 



^An error in this formula lead to some erroneous conclusuions about the j/-behaviour of the structure 
function in the original version 
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We finally come to the gluon density. Although, strictly speaking it is not a physical 
quantity, its properies have been much discussed recently in connection with its saturation 
for the large nucleus [16]. It can be related to our function (p via the standard expression for 
the structure function in its terms (see e.g [14]) 

F2ix,Q^) = / d%£r[pTir) + pLir)]Fix,r,b), (47) 

where 

and / (Pb(dxG{x, k'^,b)/dk'^) is up to a factor the gluon density in the momentum space: 

dNil) _ 1 dN{l) _ 1 f ^^^ dxG{x,l\b) ^ ^^g^ 



n dp ttJ dp 

In the following we shall study the double density in momentum and impact parameter, which 
is just dxG{x,k'^,b)/dk^. Comparing (47), (48) with the corresponding expression in terms 
of (f), which follows from (39) and (40) we find a relation 

F(x,r,b) = ^r^(i)(ln-,r,b). (50) 

TT^^ X 

Taking a Fourier transform of (48) and neglecting the term proportional to S^{k) we obtain 

dxG(x,k^,b) 2N,o„2 2n 1 , 

' ^ = — A;2vi0 In -, k, b). 51 

Ok^ TT^^ X 

This is the desired relation between our function (j) and the gluon density in the combined 
momentum and impact parameter space. 

Applying k^V\ to the found function we thus find the gluon density up to a trivial 
numerical factor evident from (51)( ~ 0.76 with = 0.2). Function h{k) = k'^\7l(j){y,k,b) 
for different values of y and B = 0.12 and 0.02 is shown in Figs. 8,9. Its form at different 
y results quite remarkable. As one observes, at any given rapidity the found density has the 
same, roughly Gaussian shape in variable ^ = InA;, centered at the point ^ = ^o(y) very near 
to ^1 at which the straight line (38) crosses the .x-axis (sec Section 2). With the growth of y 
the distribution moves to the right with a nearly constant velocity practically preserving its 
form. Approximately the distribution can be described by 

/i(A:) = /ioe""(^"^"(^))', (52) 

where ho and a are practically independent of y and ^o(y) linearly grows with it: 

ho 0.3, a -0.3 eo(?/) = eoo + 2.23y. (53) 

The only quantity which clearly depends on the initial distribution is the starting position 
^00) so that for different initial functions the picture in Figs. 8 and 9 shifts along the ^ axis 
as a whole. 

Evidently at a given value of k the density stays always limited, irrespective of the form 
of the initial distribution (and on the atomic number A, in particular). In this sense we have 
saturation as discussed in [16]. However with the growth of y the strongly peaked density 
moves away toward higher values of k so that the density at a fixed point tends to zero at 
high values of rapidity. We thus have "supersaturation" : with y — > oo the gluon density at 
an arbitrary finite momentum tends to zero. 
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Comparing Figs. 8 and 9 one can see how the memory about the initial distribution 
(except for ^oo) is gradually erased in the course of the evolution in y. For a very peripheric 
collison off lead {B = 0.02) at the initial stages of the evolution the density is correspondingly 
much smaller than for a central collision [B = 0.12). However already at y = 2 the form of 
the distribution is practically indistinguishable from the central collision, the only remaining 
difference being the shift along the ^ axis. 

5 Conclusions 

The BFKL fan diagram equation has been solved numerically in the large range of rapidities 
up to y = 50. The main results are the following. 

The idea that the fan diagrams themselves satisfy the unitarity condition has been sup- 
ported by the fact that the found cross-sections for the scattering of a qq pair off the nucleus 
tend to a constant value at high rapidities. Since the found cross-sections do not rise with 
energy, the leading j-plane singularity turns out to be a simple pole at j = 1. The limiting 
cross-sections prove to be universal: they do not depend on Q^, that is, on the transverse 
dimension of the qq pair, nor on the initial colour distribution inside the nucleon. Their de- 
pendence on the target nucleus thus reduces to a scale factor R\. Physically they correspond 
to scattering of a colour dipole off a black disk. 

The nuclear structure function does not saturate at high rapidities, due to the singularity 
of the transverse distribution prif) at r = 0, which makes it non-normalizable. At large y it 
continues to grow nearly linearly in y. 

These results fully agree with predictions made in [11] on the basis of the found per- 
turbative solution of the BFKL fan diagram equation and asymptotic estimates made in 
[12]. 

A completely novel result concerns the gluon density of the nucleus. At sufficiently high 
rapidities, greater than 10, the gluon density aquires a form of the soliton wave in y — In fc 
space, which, with the growth of y, moves along the Ink towards greater k preserving its 
nearly Gaussian shape. Thus at any finite k the gluon density eventually goes to zero at high 
enough k 

Finally we have to stress that all these properties begin to be clearly visible only at very 
high rapidities and momenta: y > 10 and k > 100 GeV/c with ag = 0.2. With smaller a's 
these values grow correspondingly. 
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7 Appendix. Colour dipole approach 

To compare our BFKL fan diagram equation to that of Kovchegov [10] we present here a short 
derivation of Eq. (1) from the colour dipole approach, which will make clear the difference 
between the two equations. 
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In the colour dipole approach the single scattering term (2) is presented in the form 

A{y,b) = -2isAT{b) J d\p{r) j d\ini{r,n,y)T{n), (54) 

where 

r{ri) = J (fr'G{0, n, r')m(r') (55) 

and ni{r,ri,y) is a single dipole density at rapidity y introduced by A.H.Mueller; ri and r 
are the dipole lengths at rapidity y and at y = respectively, so that 

ni{r,ri,y)y=o = S^{r -n). (56) 

(Note that as in [10] our ni is A.Mueller's one divided by 2-rrrf.) 

To introduce the multidipole densities A.H.Muller constructed a generating functional 
Z{ri,ro,y\u), where u = u{pi,pf) is a function of two dipole endpoints in the transverse 
space (which, for brevity, Iwe denote by a single symbol p for) . The functional Z satisfies the 
following nonlinear equation 

Z{ri,r2,y\u) = u{ri,ro)e^y'^^'^°^ + 

^2 AT ry p J.2 

^ / d\2^^Z{n,r2,y'\u)Z{r2,ro,y'\u), (57) 

where rio = ri — vq etc and Lo{r) is the gluon trajectory u;(q) in which the momentum q is 
substituted by r: 

^W = -^ln- (58) 

with e the cutoff at small r (in the ultraviolet). The functional Z is normalized according to 

Z{ri,ro,y\u)u=i = l. (59) 

The A;-fold inclusive dipole density is given by a fc-fold derivative of D with respect to u at 
u=l: 

1 Z 

nk{ri,ro,y;pi,...pk) = tt c / x — . , . ■ (60) 

k\ 6u{pi)....5u{pk)u=i 

At this point we make our first comment as to the comparison with [10]. Our form of the 
functional equation for Z is the same as in the original pater of A.H.Mueller [1] and in [10], 
except for a slightly different choice of dipole coordinates and for the order of arguments in 
the 2nd Z in the non-linear term: their form would correspond to Z{rQ,r2,y'\u). Comparing 
with the BFKL equation for the single dipole density, one can verify that our choice is better. 
However this point is irrelevant for the following, since in the interaction with the nucleus 
only even functions of rio appear. 

The dipoles are to interact with the nucleus target with a zero transferred momentum. If 
the dimension of the dipole is smaller than the internucleon distance in the nucleus then it 
will interact with a single nucleon as a whole. This picture lies at the basis of the standard 
fan diagram approach, where each pomeron finally interacts with a single nucleon. We used 
precisely this picture in our derivation of Eq. (26) in Section 2. In this picture the only trace 
of the nucleus will be an additional factor AT(b) multiplying the interaction with the nucleon 
r(r) (Eq. (54)) at a given impact parameter b. In [10] a different idea is exploited: it is 
assumed that each of the dipoles can interact with many nucleons. The latter interaction is 
assumed to have an eikonal form. So the interaction with the nucleus they consider is a two- 
stage one: first the projectile generates many colour dipoles (BFKL fan diagrams) and then 
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each dipole multiply interacts with the nucleus a-la Glauber. Although technically it is not 
difficult to take into account the final eikonalization of the interaction, appropriately changing 
the single scattering term r(r) in (55) and the following formulas, we do not think it is 
reasonable. On the one hand, as we shall sec, in the interaction with the nucleus the densities 
are considerably damped at large distances as compared to the usual BFKL behaviour. The 
confinement should further restrict their spatial dimensions. So for a nucleus with a large 
internucleon distance it does not seem reasonable to assume simultaneous interaction of a 
dipole with two or more micleons. On the other hand, should such interactions be really 
important, one cannot expect to correctly describe the interaction with the nucleus of a 
dipole of a given (and fixed) dimension by the Glauber formula. Note that the expression 
used in [9] for it does not correspond to the BFKL picture for the scattering on a single 
nucleon. So we take r as given by (55) and this is the main difference between our derivation 
and that of [10]. 

With a chosen r, the interaction with the nucleus will be described by densities 

i^k{ri,ro,y) = j fl {d'^PjT{Pj)AT{b)jnk{ri,ro,y; pi, ...pk)- (61) 

Here r(p) = r(p/ — pi) depends only on the dipole length and is given by (55). 
Differentiating (57) and using (61) one easily obtains 



- I ay e ' ' j a r2- 

''l2'20 



^ 3 ^dy'e2^no){J/-J/') / dV24V[^^i(n,r2,y') + i^i(^2,ro,2/')] (62) 



and for > 1 



Uk{ri,ro,y) = ^ / / VVK(^1' ^2, y') + ^^fe(^2, ro, y') 

Stt-^ Jo j rinrnn 



12' 20 

k-i 



+ Yl ^ji'^i^r2,y')j^k-j{r2,ro,y')]- (63) 

i=i 

We suppress the evident dependence on the impact parameter b. From the structure of the 
equations and the form of the inhomogeneous term it follows that the densities i^(ri,ro,y) 
depend only on the initial dipole length no. As a result, the two terms separated from the 
sum over j on the right-hand side give the same contribution. One then finally finds equations 

^2 Af ry f ^2 

Mno, y) = ^r(fe)r(rio)e22'-(^io) + ^ / dy'e2-('-io)(^-^ ) / dVs^V ^1(^20, y') (64) 



and for A; > 1 



-.(no,y) = ^ rdi/'e2-(-o)(.-.') [ d\,^Mr20,y')+ 

47r Jo J '^\2'^20 



^ r dy'eMr.o)iy-y') j g -Ari2,y>k-Mo. v')- (65) 

•> "12^^20 j=i 

The total forward scattering amplitude on the nucleus will be given by the expression 
(54) in which the single dipole interaction / d^rini(r, ri, y)r(ri) is substituted by the sum of 
all multidipole interactions '}2k ^k{^)- Presenting the amplitude in the form (10) we have 

^{r,y) = -Y.Uk{r,y). (66) 

k=l 
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AT ny n ™2 



Summing Eqs. (64) and (65) over k we obtain an equation for 

1 / " W-,/^2a;(^io)(?/-?/') f T 

47r3 

^ r d2;'e2-('-")fe-^') / dV24iV^(n2, y'Mr^o, y'). (67) 

Comparing this equation with the one derived in [10], apart from a different inhomoge- 
neous term, which was discussed earlier, we find difference in the spatial arguments of the 
function Our $ depends only on one such argument: the dipole dimension ri2. The 
equivalent Kovhegov's function depends on two spatial arguments: it depends not only on 
the dipole dimension but also on its center-of-mass coordinate bo, not to be confused with the 
impact parameter b which does not enter his equation at all. The dependence on bo should 
be governed by the inhomogeneous term, which seems to be independent of bo (Eq. (6a) of 
[10]). Then the dependence of on the 2nd argument seems to disappear and Kovchegov's 
equation coincides with (67). However this contradicts his initial formula (Eq. (4) of [10]) 
in which one integrates over all bo- Modulo all these (small) inconsistencies and a different 
inhomogeneous term, our Eq. (67) coincides with Kovchegov's. 
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9 Figure captions 

Fig. 1. The equation summing fan diagrams. Lines represent pomerons. 

Fig. 2. (j) as a function of momentum at different y (in units n/asN) for central collisions on 

lead {B = 0.12, Eq.(359)). 

Fig. 3. (j) as a function of momentum at different y (in units n/asN) for peripheral collisions 
on lead {B = 0.02, Eq.(35)). 

Fig. 4. Cross-section (with — > 1) as a function of y (in units tt/usN) for different Q^. 
Fig. 5. Cross-section aL (with — > 1) as a function of for different y (in units ir/agN). 
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Fig. 6. The structure function of lead F2a{x,Q'^) as a function of x for different Q^. 
Fig. 7. The structure function of lead F2a{x, Q^) as a function of for different x. 
Fig. 8. The gluon density (in units 7r'^N/2as) as a function of momentum at different 
rapidities y for central collisions on lead {B = 0.12, Eq. (35)). 

Fig. 9. The gluon density (in units 7r'^N/2as) as a function of momentum at different 
rapidities y for peripheral collisions on lead {B = 0.02, Eq. (35)). 
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